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The Hubbard model on the honeycomb lattice describes charge carriers in graphene 
with short range interactions. While the interaction modifies several physical quan- 
tities, like the value of the Fermi velocity or the wave function renormalization, the 
a.c. conductivity has a universal value independent of the microscopic details of the 
model: there are no interaction corrections, provided that the interaction is weak 
enough and that the system is at half filling. We give a rigorous proof of this fact, 
based on exact Ward Identities and on constructive Renormalization Group methods. 



1. INTRODUCTION AND MAIN RESULTS 



The effects of interactions in quantum many body theory at low temperatures pose noto- 
riously difficult problems; in certain cases the physical properties of the system are radically 
changed with respect to the non interacting situation, while in other cases, like in the so- 
called Fermi liquids, a simple modification or renormalization of the physical quantities is 
expected. There is, however, a very small group of phenomena which are universal; the 
physical quantities appear to be protected from any renormalization due to the interactions 
and their values do not dependent on the details of the model, but rather upon fundamental 
constants. A celebrated example is provided by the Quantum Hall Effect, in which the value 
of the plateaus only depend on the von Klitzing constant h/e'^ and not on the material pa- 
rameters; universality appears to be related to topological invariance [1, 2] or to the presence 
of Ward Identities [3]. Other examples of universal phenomena come from the physics of 
superconductivity, in which the magnetic quantum h/e plays an important role. 
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In recent times, evidence for universality has been observed in the conductivity of 
graphene, a one atom thick layer of graphite. Its electronic properties can be well described 
in terms of a tight-binding model of electrons hopping from one site to a neighboring one 
of a honeycomb lattice, but often this model is approximated by an effective one expressed 
in terms of massless Dirac fermions in the two-dimensional continuum [4] . Recent optical 
measurements [5] show that at half-filling and small temperatures, if the frequency is in a 
range well inside the temperature and the band-width, the a.c. conductivity is essentially 
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constant and equal, up to a few percent, to ctq = y|. Such value only depends on the 
von Klitzing constant and not on the material parameters, like the Fermi velocity; it is 

apparently universal, at least inside the experimental precision. 

Is graphene a.c. conductivity truly universal? Theoretically, the computation of the 
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conductivity in the absence of interactions gives exactly the value ctq — both in the 
idealized Dirac description [6] and even in the more realistic tight binding model [7]. How- 
ever, since truly universal phenomena are quite rare in condensed matter, it is important 
to understand whether this apparently universal value is just an artifact of the idealized 
description in terms of non-interacting fermions or rather it is a robust property still valid 
in the presence of electron-electron interactions, which are certainly present and expected 
to play a role in real graphene. Such a question has been studied in the physical litera- 
ture, but contradictory results have been found, see [8-10]. The reason for this is that the 
Dirac approximation, which works well for the free gas, is not very accurate in the pres- 
ence of interactions: the corrections to the conductivity are expressed by logarithmically 
divergent integrals. One can argue that such divergence is spurious, just an artifact of the 
Dirac approximation, and that a regularization must be adopted to cure it (in the tight 
binding model, the lattice provides a natural cut-off); however, the results appear to be 
regularization-dependent and no unique predictions can be drawn. 

In this paper we consider the Hubbard model on the honeycomb lattice, as a model of 
monolayer graphene with screened interactions, and we prove that the a.c. conductivity has 
the universal value o"o = xf ^^^^ presence of interactions: the interaction corrections 
to the conductivity are vanishing, provided that the interaction is weak enough and the 
system is in the half filled band case. Remarkably, the presence of the lattice and its 
symmetries is essential to get the result. The idea of the proof is based on the two main 
ingredients: (i) exact lattice Ward Identities (WI) relating the current-current, vertex and 
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2-point functions; (ii) the fact that the interaction-dependent corrections to the Fourier 
transform of the current-current correlations are differentiable with continuous derivative 
(in contrast, the free part is continuous and not differentiable at zero frequency). 

The paper is organized in the following way. In Section 1 A wc describe the model, in 
Section 1 B wc derive the Ward Identities, in Section 1 C (and Appendix A) we perform 
the computation of the conductivity in the non interacting case, and in Section 1 D we 
present the proof of our main result, under some regularity assumptions on the current- 
current correlations. The proof of these regularity properties, the full description of the 
expansion for the current-current correlations and vertex functions, as well as the proof of 
convergence of this expansion, is given in Section 2 (Appendix B and C collect the proof of 
some symmetry properties extensively used in Section 2). 



We consider electrons on a two-dimensional honeycomb lattice interacting via a local 
Hubbard interaction, as a model describing the charge carriers in graphene. Its ground 
state properties at half-filling, including the asymptotic behavior of the correlations at 

large distances, have already been analyzed in [11, 12]. In this subsection we recall the 
definition of the model and introduce some of the key observablcs (density and current), 
which will allow us to define the conductivity, i.e., the quantity of main interest in this paper. 

— * — * 

The fermionic fields and the Hamiltonian. Let A = {nili+n2l2 '■ ni,n2 — 0, ■ ■ ■ , L—1} be 
a periodic triangular lattice of period L, with basis vectors: h = |(3, \/3), I2 = |(3, — \/3). 
Let us denote by A^ = A and A^ = A + 5, the A- and i?-sublattices of the honeycomb 
lattice, with 5i the n.n. vectors defined as: 



If A and B are two arbitrary constant vectors, we introduce creation and annihilation 
fermionic operators for electrons sitting at the sites of the A- and B- sublattices with spin 



A. 



The model and the observables 




(1.1) 
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index a as 

4,. - E ' ^^^^^ (1-2) 

where 5^ = = mGi/ L + n2G2/ L : < < L}, with Gi,2 = f (1,±V3), is the 
first Brillouin zone; note that in the thermodynamic hmit -^"^ X^ggg^ — > \B\~^ j^dk, with 
B = = ^iGi + : 6 e [0, 1)} and \B\ = 87rV(3v^). The operators ai^, fei^ satisfy 
the canonical anticommutation rules and are periodic over A; their Fourier transforms are 
normalized in such a way that, ii k,k' are both in the first Brillouin zone: 

4, J = L\^A,-e'Sa,.' , b< J = L\-^A,-e'Sa,a' ■ (1-4) 

Moreover, they are quasi-periodic over the first Brillouin zone: 



k+Gi,a k,T k+Gi,T k,T 



where e^'^*^^ — e*^"^/^, for all values oii,j. The phases A and B are arbitrary, the freedom in 
their choice corresponding to the freedom in the choice of the origins of the two sublattices 
and A^ (this symmetry is sometimes referred to as Berry-gauge invariance). A convenient 

choice for A and B, which makes the fields a- , hi periodic over the reciprocal lattice A*, 

k,a k,a 

is ^4 = and B = 6i, which reads: 



4,a 



^-'E-""«t' = E ' -e^- (1-6) 

This is the choice made in [11, 12], which will be used throughout this paper, too. 

The grand-canonical Hamiltonian of the two-dimensional Hubbard model on the hon- 
eycomb lattice at half-filling is H\ = Hl{t) + UV\, where H^{t) is the free Hamiltonian, 
describing nearest neighbor hopping e M is the hopping parameter): 

Hlit) = Hint,,}) , Hint,,}) = - E E + %&t^,^,.%,.) , 

j=l,2,3 

(1.7) 

and V\ is the local Hubbard interaction: 

= E (4t%,t - 1) - ^) + E (^tt^5,t - 1) i^tA - 1) ■ (1-8) 

1=1,2,3 1=1,2,3 
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The current and density operators. The current is defined as usual via the Peierls sub- 

— * 

stitution, by modifying the hopping parameter along the bond {x, x + Sj) as 

t tgj{A) = t e'^-^o Ax+sS^ySj ds ^ (^_g) 

where the constant e appearing at exponent is the electric charge and A{x) e is a 

— * — * 

periodic field on S\ — {x — L^ih + L^2h '■ £ [0, 1)}. Its Fourier transform is defined as 
A{x) = \Sa\-' J2pev^ Ape-'^, where \Sa\ = and = {p = mGi/L + riaGa/L : m e 

Z}; note that in the thermodynamic limit |iSa|~^ 12peVA ~^ Jr^ dp. If we denote by 

Ha{A) = H^{{t£j{A)}) + UVa the modified Hamiltonian with the new hopping parameters, 
the lattice current is defined as J^^'' — —\S\ \ dH\{A)/dAp, which gives, at first order in A, 



geVA 



where, if rji ■ 



P 1,2 Z-^ •? 'pV k+p,a k,<7 k+p,a k,a ^ 

keBL 

o-J i=l,2,3 

is the paramagnetic current (in the last rewriting, vo — St/ 2 is the free Fermi velocity and 
Ji — (2i/3)(at bZ ?• — bt ■? aZ) are the bond currents) and 

X \ I / \ x,a x+dj,a x+Sj,a •'^i"^ ' 

[^P,^]lm- E e-^^^^^^'{SMS,)rnriivi^x,j , (1.12) 

xeA 

j=l,2,3 

with A^,- = — e^tV (at h" ^ + 6^ ^ ), is the diamaqnetic tensor. Similarly, the 
density operator is defined by coupling the local density to an external field, i.e., by adding 
to the Hamiltonian a local chemical potential term of the form 

ma{h) = - E E /^(^)4,.«.> -EE ^(^HaK^a , (1-13) 

where ii{x) — |»5a|~^ ^-^^^ //^e"*^^ is a periodic field on Sa. If we denote by i?A()u) = 
H^{t) + UVa + Ma{ijl) the modified Hamiltonian in the presence of the local chemical po- 
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tential, the lattice density is defined as p^— —\S\ \ dH\{p)/dpL^, which gives 

= J- y (at ^ + e-'^H± . hz ) ^ (1.14) 

a=n 

X^Aa xeAg X^Aa xeAg 

It will be convenient for the incoming discussion to think the two components of the param- 
agnetic current J^^i, I — 1,2, as the spatial components of a "space-time" three-components 
vector J^^n, /i — 0,1,2, with J^^ — epp. In the following, it will also convenient to introduce 
the reduced current j^, related to the paramagnetic current by 

Jp^VQjp, (1.15) 

with vq = 3t/2 the free Fermi velocity. 

Schwinger functions and response functions. The thermal state of the system at inverse 
temperature (3 > 0, associated to the density matrix e~^^'^, can be characterized in terms 
of Schwinger functions and response functions. The Schwinger functions are the analytic 
continuation to imaginary time of the off-diagonal elements of the reduced density matrices; 
the response functions give us informations about the reaction of the system to a diversity 
of external probes within the linear response regime. They are defined as follows. Let 
ib^^ — (at , bt ^ ), let i = 1, . . . , n, be local monomials in the ib^^ operators and let 
US denote by Ox^ = QXi,oHAQ(^)Q-xi,oHA ^j^g corresponding imaginary-time evolved operators; 
here Xj = {xi^, Xi) and xi^ G [0, (3) is the imaginary time. The average of a product of local 
operators in the thermal state of the system at inverse temperature /? > is defined as 

,Oa)...o<»)> ^ Tt{e-^"-'T(0g...0g)} 

where T is the operator of fermionic time ordering, acting on a product of fermionic fields 
as: 

where Si e {+, -}, Gi e {t,4}, Pi e {1,2}, e [0,^) X A, Vx,<,,i = <^ and, if Si = (0,5^, 
V'x,o-,2 — ^x+Sicr- Moreover, tt is a permutation of {l,...,n}, chosen in such a way that 
3^7r(i)o > • • • > 2:^(n)o, ^ud (—1)'^ is its sign. [If some of the time coordinates are equal each 
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other, the arbitrariness of the definition is solved by ordering each set of operators with the 
same time coordinate so that creation operators precede the annihilation operators.] Finally, 
we denote by (O^i', ■ ■ ■ ', O^}) ^ ^ the corresponding truncated expectations. We shall also use 
the notation {■) ^ = limL_,oo and (■) = limi^oo (■)/3- 

Choosing the local operators O^j in Eq.(1.16) simply as monomials in the fermionic fields, 
we get the Schwinger functions of order n: 

Choosing the operators 0^} as suitable combinations of the current and density operators, 
we get the current-current, density-density and current-density response functions: 

, 1-13/2 .13/2 

^M^(p) -Rn / dyo^-''"'^^°-'''\JM,,; Jiyo,-p),^)pL ' (1-19) 

pLi J-13/2 J-13/2 

where p = {po,p) G x Va and J{xo,p),fi = e^"^'^ Jp,^e~^°^^. An important role will be 
also played in the following by the two- and three-points functions: 

^'''^(k) 4i / d^f dye^''^^--y^S^'\^,-,a;y,+,a), (1.20) 

J{/3,L) J (0,1) 
1 p(3/2 1-0/2 .0/2 

(1.21) 

where /(^,^)C?x is a shorthand for fJ^j^dx^Y^i^K and = e^o^AV'fe,^^"''"^^, with 

■0^ = (ci^ , 6^ ). Here and in the following, we will exploit whenever possible the vectorial 
structure of and the tensorial structure of products of V^'s; to this purpose, we will think 
of as a column vector and ip^ as a row vector, so that, e.g., ip^„ipy„ will be naturally 
thought as a scalar, while V'xo-V'y o- will be thought as a 2 x 2 matrix. In particular, both 
^'^'^(k) and G2,iJ^(k, p) (for any fixed choice of e {0, 1, 2}) can be thought as 2x2 matrices. 

Conductivity. The ac conductivity in units such that ^ = 1 is related to the current-current 
correlations via Kubo formula [7], which reads, for all 27r^~^Z 9 po 7^ and Z, m e {1, 2}: 

'^L(Po) = [KLiPo, 0) + AfJ , (1.22) 

where i^L(P) = limL^+oo i^f^ (p) and 

= iToo 1? ^ (^.^^^^-(A.-,),,^ (1.23) 
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is the diamagnetic contribution (see Eq.(1.12)) and the factor 2/(3-\/3) appearing in Eq.(1.22) 
must be understood as the inverse of the area of the unit cell of the hexagonal lattice. 

The main goal of this paper is to compute C;^j(po) in the zero temperature and zero 
frequency limit (taking the limits in a suitable order, so to make contact with experiments 
on the optical conductivity of graphene), i.e., to compute the so-called universal optical 
conductivity: 

aim lim hm (7f^(po) ■ (1.24) 

A key role in its computation will be played by Ward Identities, which show that the 
quantities introduced above arc not independent; on the contrary, they are related by exact 
identities, which we now describe. 

B. Conservation laws and Ward Identities 

By definition of p{xo,p) = e^°^^p^~^°^^, we have: 

dxoP{xo,p) = [Ha,P{xo,p)] ■ (1-25) 

Computing explicitly the r.h.s. of this equation and using, in particular, the fact that 
[Va, Pf] = 0, we find: 

d.oPixo,P) = t e-'^^il - e-'^'^)ial,bZ^s„a - bt+s„a<^Z,a) ■ (1-26) 

Comparing the r.h.s. of this equation with the definition of the paramagnetic current, we 
recognize that Eq.(1.26) can be rewritten as a continuity equation: 

-iedxoP{xo,p) + ip- J{xo,p) = . (1.27) 

Using the continuity equation Eq.(1.27), we can easily derive an exact identity relating 
three- and two-point functions. In fact, by the definition of G2,'i^Q(k, p), see Eq.(1.21), and 
integrating by parts, we find: 

ri3/2 .13/2 

(1.28) 
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The derivative with respect to Zq can act either on J(^zo,p),o-, in which case we use the continuity 
equation, or on the Heaviside step functions involved in the definition of time-ordering 
operator T. After some straightforward algebra, and using the fact that 

= -"^m,. ' - n^p,.^^ ' (1-29) 

with 




(1.30) 



we end up with 



-ipoG&(k,p) +PiG5i':i(k,p) +p2GC'i':2(k,P) = -e^^'^(k + p)M(p-) + eM(p-)^^'^(k) . 

(1.31) 

Proceeding in the same way, we also find: 



-ipo^o^;o^(p) + Pi^f;o^(p) + P2^|/ (p) = , (1.32) 
-ipoi^o^;^ (p) + Pii^f,;^(p) + P2i^2^i (p) = -T^ ■ {K^-p)p,i\ , (1-33) 

where m — 1,2, and we used the fact that 

e[Pp, J-p] = (p • ^j) Sj\rii\'^Agj = p- Ap-_^ . (1.34) 

The term in the r.h.s. of Eq.(1.33) is known as Schwinger term. 

Note that from Eq.(1.33), setting, e.g., p2 = 0, we find that, for i = 1, 2, 

^h(Po,Pi,0) + ^hm^-^([A(pj,o),(-pi,o)]iJ^^^ = i^i^oi(Po,Pi, 0) , (1.35) 

provided the limits of these functions as ^, L — > oo exist. 

Now, if we knew that the correlations in Eq.(1.35) were continuous and continuosly differ- 
entiable p = 0, from the latter equation we would conclude that limp^o[-^ij(p) + — 
(with = lini/3-s.oo sec Eq.(1.23)); moreover, continuous differentiability together 
with the symmetry properties of our model would imply that (with a proof analogous to the 
one of item (ii) in Proposition 1 below) limp^o f^po-^ij(p) — 0- Comparing these equations 
with the definition of the universal optical conductivity, see Eqs.(1.22)-(1.24), then we would 
be tempted to conclude that Uij = (and a similar argument would imply that (72i — 0). 
This is in contrast with the explicit computation of the conductivity in the non-interacting 
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case, which will be discussed below. The solution to this apparent paradox is that Kim{p) is 
not continuosly differentiable at p = 0. In fact, as it will turn out from the following discus- 
sion, the regularity properties of the Fourier transform of the current-current correlations 
play a crucial role in the physical properties of the conductivity 



C. Properties of the non interacting model 

In the absence of interactions, that is for U — 0, the two-points function defined in 
Eq.(1.20) reads (see [11, 12]): 

-S^'^ k) = ^ \ ° ^ M =: So{k) , 1.36 

where vq = |t and Q{k) = | Y^^=i^^'^'"^^~^^^ ■ The complex dispersion relation Q{k) vanishes 
only at the two Fermi points 

P/ = (|.±^). (1.37) 

close to which it behaves as follows: ^{p^ + k') = ik[±k2 + 0{\k'\'^). The Schwinger functions 
of higher order can be explictly computed in terms of 5*0 (k) via the fcrmionic Wick rule. 

Also the universal optical conductivity can be computed explicitly, see Appendix A, and 
turns out to be equal to (restoring the presence of the dimensional constant H — /i/(27r) in 
the result): 

I 2 2eX y fdko f dk^ ^ 5o(k+(po,0))-5o(k) ^,^^,^ ^^>^ 

"^^■1-0 = ^^.o^"^J ^ LW\ i Vo r.(A;,0)5o(k)r,(A;,0)| , 

(1.38) 

where 

^ / ie-ikSj-Si)\ 

Note that the r.h.s. of Eq.(1.38) depends on vq. Moreover, the integral is not uniformly 
convergent in po as Pq — )■ 0; in particular, it is well known that one cannot exchange the 
limit with the integral [13]. An explicit computation, see Appendix A, yields 

a value that, remarkably, does not depend on vq- It is also remarkable that the same value 
of the conductivity is found in the so-called Dirac approximation, that is by replacing 5'o(k) 
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in Eq.(1.40) with its linear approximation around the Fermi points, in the presence of an 
ultraviolet cutoff, i.e., 

x-^ ^ „ ^ - / —iko —vo(—iki+u!k2) 

uj=± \-VQ{iki + ujk2) -iko 

(1.41) 

where £ is a small positive number. Therefore, (Jij\^_^ does not depend on the lattice 
parameters and, in this sense, we can say that the free conductivity is universal. 

D. Universality of the conductivity 

Let us now discuss the effects of the Hubbard interaction. We will first recall the results of 
[11, 12], concerning the thermodynamic functions of the model, the Schwinger functions and, 
in particular, the two- and three-points functions. Next, we will state our results concerning 
the interacting conductivity. 

In Theorem 1 of [11], we proved that, if U is small enough (uniformly in the system 
size and in the temperature), then the specific free energy fii{U) and the finite temperature 
Schwinger functions in the thermodynamic limit are analytic functions of [/, uniformly in 
/3 as ^ — >■ oo, and so are the specific ground state energy e{U) and the zero temperature 
Schwinger functions in the thermodynamic limit. In particular, we proved that the Fourier 
transform of the zero-temperature two-points function in the thermodynamic limit is singular 
only at the Fermi points k = = (0,p^) and, close to the singularities, if a; = ±, it can 
be written as 

1 I -iko -VFn*(p'^ + k')\ ( ,\ 
5(p^ + k') = - ° . M (l + i?(k')), (1.42) 

where Z and vp are two analytic functions of C/, analytically close to their unperturbed 
values, 

Z = l + 0([/'), vf^Vo^O{U'^) . (1.43) 

Moreover the matrix i?(k') satisfies ||i?(k')|| < C|k'|'' for some constants C, > and |k'| 
small enough. In [12] we also announced the following result, whose detailed proof will be 
given below. 
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Theorem 2. There exists Uq > such that, the three-point function Eq.(1.21) is uni- 
formly analytic in j3,L and uniformly convergent to an analytic function as (3, L ^ oo, for 
aWk ^ Pp, p and for (min^ |k — p'^|), |p|-(min(^ |k— p'^l)^-'^ sufficiently small. Moreover, 
the vertex function at the Fermi points (in the thermodynamic and zero-temperature limits) 
is equal to: 

liin lim [^(k + p)]-'G2,vA^, p)[S{k)]-' = eZ.r.if^p, 0) , (1.44) 

where Vi{k,p), i = 1,2, were defined in Eq.(1.39), To{k,p) = 1 and = Zn{lJ), /x = 0, 1, 2, 
are analytic functions of U . 



Note that the existence of the hmits in Eq.(1.44) is part of the statement of the theorem. 
An important consequence of Theorem 2 (also announced in [12]) is obtained by combining 
its result with the WI Eq.(1.31): in fact, using Eqs.(1.42)-(1.44) into Eq.(1.31) we find that 
the vertex functions are related to the wave function renormalization Z and to the Fermi 
velocity by simple identities: 

Zo = Z , Zi = Z2 = vpZ . (1.45) 

These relations can be proven as follows: take the limits in the l.h.s. of Eq.(1.44) with 

k = p'p + k'e^ and p = p'e^, where Gq = (1,0,0), Si = (0,1,0) and 82 = (0,0,1); using 
Eq.(1.31) and Eq.(1.42) we can rewrite the l.h.s. of Eq.(1.44) (after having taken p' 0) 
as: 

-e(^)'^-" lim[^(p:^ + k'e^)]-^d^S{p^p + A;'e^)[^(p:^ + k'e^)]-^ = 
«'— >^o 

= e(i)^-° ^^d,[S{p^p + k'e^)\-^ (1.46) 

k'— >^0 

that, using again Eq.(1.42), is equal to eZ{vFY~^^'''T 



To summarize, at half filling, weak local electron-electron interactions do not change the 
infrared behavior of correlations: they "just" change the values of some physical parameters, 
namely the wave function renormalization Z, the Fermi velocity vp and the vertex functions 
Z^; the latter are related to Z and vp in a simple way, thanks to Wis. The next natural 
question we would like to answer to is how is the conductivity changed by the presence 
of interactions. Since the infrared behavior of the interacting correlations is the same as 
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the non-interacting one, it is natural to expect that the interacting conductivity remains 
finite in the zero-temperature and zero- frequency hmit; what is apriori unclear is whether 
the zero-frequency conductivity remains universal in any reasonable sense, in analogy with 
the universal behavior of the free conductivity. Quite remarkably, wc can prove that the 
interacting conductivity is universal in a very strong sense: namely, we prove that aim is 
not only independent of the details of the lattice, but it is also exactly independent of U. 

Theorem 3. There exists a constant Uq > such that, for \U\ < Uq and any fixed po 

(nan vanishing and sufficiently small), (JimiPo) is analytic in U uniformly in j3 as /3 —)■ oo 
and uniformly convergent to an analytic function of U as (3 ^ oo. Moreover, 

aim = lim lim af^{po) = -rT^Sim ■ (1-47) 

Note that the limit /3 — >■ oo is taken before the limit po — >■ 0"*". In other words, the 
theorem says that the interaction corrections to the conductivity are negligible at fre- 
quencies (3"^ <^ Po <^ t, in agreement with experiments on the optical conductivity [5]. 
The above results says that all the interaction corrections to the conductivity cancel out 
exactly, even if the Fermi velocity and the wave function are renormalized by the interaction. 

The proof of Theorem 3 is based on two main ingredients: (i) the use of the exact Ward 
Identities Eqs.(1.32)-(1.33); (n) the fact that the interaction-dependent corrections to the 
Fourier transform of the current-current correlations are differentiable with continuous 
derivative (in contrast, the free part is continuous and not differentiable at zero frequency). 
The main technical point of this paper is to control the regularity properties of the in- 
teraction corrections to the conductivity, which are summarized in the following proposition. 

Proposition 1. There exists Uq > such that, if \U\ < Uq, then the current- current 
function Ki^{p) is analytic in U, uniformly in (3,L, for all sufficiently small p 7^ and 
it is uniformly convergent as (3, L to the function Kim{p) — lini/3_).oo liniL-)-oo ^imip)) which 
is also analytic in U for all sufficiently p 7^ 0. The latter function satisfies the following 

properties: 

1. Kim{v>) is continuous for all sufficiently small p e M x B (in particular atp — 0) and 
continuously differentiable for all sufficiently small p 7^ 0. 
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2. Kim{p) can be decomposed as: 



(1.48) 



where: 



(i) if the average with respect to the density matrix e ^^a^I'^p') associated to 

the non-interacting Hamiltonian with Fermi velocity vp, then 



-, 1-13/2 1-13/2 
P,L-^oo PL 7-/3/2 J-^/2 



L ' 



with j(^xo,p),i = G^°^'^jp,iG and jp the reduced current defined in Eq.(1.15); 

(a) R{v) is continuously differentiable for all sufficiently small p e M x B (in partic- 
ular at p — 0); moreover, R{po, 0) = R{—po, 0). 

One can immediately realize that Theorem 3 is a simple corollary of Proposition 1 and 
of the Wis Eqs.(1.32)-(1.33). In fact, from Eq.(1.33) computed at p = (0,pi,0), we find 



Ku{0:Pl: 0) + ^hm^ j2([Kim-PuO)] u) 0,L = ^ 

that implies, using the continuity of Kim{p) at p = stated in Proposition 1, 



:i.49) 



lim Kii{p) 

p->0 



:i.5o) 



and a similar argument shows that 



lim Ki^{p) = - Aim 
p— >o 



(1.51) 



for all l,m e {1,2}. Therefore, using again the continuity at p = of the current-current 
function and the definition of conductivity, we can rewrite 



lim 



3V3 po^o+ po 

We now use the decomposition Eq.(1.48) to rewrite the latter equation as 



(1.52) 



lim 



1 ZlZn 



^•^(po,0),/'-^-(po,0),m) ' ^ {joXiJO, 



\0,vp 



lim 



SVS po-i-o+ Po 



RlmiPO,0) - Rlm{0) 



;i.53) 
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Now, using the identity — Zvp, see Eq.(1.45), we conclude that the hmit in the first 
fine reduces to the computation of the free conductivity Eq.(1.38) with vq replaced by vp: 
however, since the result does not depend on the Fermi velocity, from the first line we simply 

get {e^ /h){'K /2)5im. On the other hand, the limit in the second line, using the continuous 
differentiability of i?(p), reduces to — ^^(?oi?(0), which is zero, simply because i?(p) is even 
in p. 

This concludes the proof of Theorem 3, once that Proposition 1 is given. The rest of the 
paper will be devoted to proofs of Theorem 2 and Proposition 1. 

2. REGULARITY OF THE CURRENT-CURRENT CORRELATIONS 

In this section we prove Theorem 2 and Proposition 1. We will use an extension of the 
method discussed in [11]. We will assume the reader familiar with the proof in [11] and 
we will only describe in detail the new aspects of the construction, as compared to the one 
in [11]. Still, we will try to be as self-consistent as possible, possibly giving reference to a 
well-defined section of [11] for the few technical aspects that will not be fully reproduced 
here. 

A. Grassman Integral representation for the correlation functions 

The goal is to prove analyticity in U and regularity in k, p of the three-point and current- 
current functions. We remind the reader that the proof of analyticity and the control of 
regularity of the Schwinger functions (in particular, of the two-points function) has already 
appeared in [11], see in particular Section 3.4 of [11]. The starting point of our construction 
is a representation of the generating function for correlations in terms of a Grassmann 
functional integral, completely analogous to the one used in [11] to write the generating 
function for the Schwinger functions. The Grassmann functional integral we are intersted 
in is defined as follows. 

Let M e N and Xo(^) ^ smooth compact support function that is 1 for i < ao and 
for i > 2ao, with oq a constant that can be chosen equal to, e.g., 1/3 (see the con- 
dition on ao appearing after Eq.(3.41) of [11] and read it for 7 = 2). Let = 
[2ttP~^{Z + I) n {ko : xo(2~^|/co|) > 0}} x Bl. We consider the finite Grassmann algebra 
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generated by the Grassmannian variables {^k,o-,p}kes*'/~^'^ ^ Grassmann integration 
/ [ IlkeB^ ^ nCT=tf ^^k,CT,p'^^k,cr,p] defined as the hnear operator on the Grassmann algebra 
such that, given a monomial ^"'") in the variables ^kcr^, its action on Q{'^^, ^"'") is 

except in the case = HkeB* 11^=14 ^ko-p^ko-p' ^'^ ^ permutation of the 

variables. In this case the value of the integral is determined, by using the anticommuting 
properties of the variables, by the condition 

p=l,2 p=l,2 

n iidH,.jK^,p] n n = 1 (2.1) 

Let us define the free propagator matrix as 

S.-M2-m( -"""'(^'V' (2.2) 

y-t'ofi(A;) -iko J 

and the "Gaussian integration" P{d'^) as 

P{d^) = \ Tl "^'^'^^°^^"'!'^°'^^' ^^,^.l^^k.l^^,^.2^^k.2 



exp { - (/3L2)-^ 5^ ^£.^k'^k,.}- (2.3) 



Let us also introduce the generating function 

>Vm,al(A, 0) = log J P(rfvI/)e^W+(*'<^)+(A^) (2.4) 



where 



V{^)^-UJ2[ d^Kt.K.A.'^Z,,,. ^ (2-5) 
<t>) = Y. I (*x,-'^x,. + 0X.*x,.) , (2.6) 

(^, J) = y2 [ d^A^^rJ: +V0 T [ d^A^jji , (2.7) 

r=±'^(^'^) ,=1,2,3 •^(^.^) 

and, in the last line, denoting by (Xi, (T2, 173 the standard Pauli matrices, 

'■^to)' ^'^t'^' ^-C-i) 
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and defining n± — {1 ± (73)/2 and a± — (cti ± i(T2)/2, 

2i 



3 

(T=t-j. (J=U 

We will be particularly concerned with the three-points functions: 

G5f>;x,y) = e^^^^>v„,„(A^)l^^^^„, «e{+,-.1.2,3}, (2.10) 

and with the current-current response functions: 

Zy^(x;y) = e^ >Vm,/^,l(A0) ^ ^ , «, b e {+, -, 1, 2, 3} . (2.11) 

The connection between these functions and the corresponding objects evaluated in he 
Hamiltonian model of Section 1 A is provided by the following proposition, which is the 
analogue of Proposition 1 of [11]. 

Proposition 2. For any P,L < +oo, assume that there exists Uq independent of (3 and 
L such that the three-points functions G^f^'^{z; x, y) and current- current response functions 
^^^^'^'■^(x; y) at distinct space-time points are analytic in the complex domain \U\ < Uq, 
uniformly convergent as M ^ oo. Then, if \U\ < Uq and z,x, y are three distinct space- 



time points, 



e(Pz;V'x,aV'y,a)^,i = Ji^G'?if/(z;x,y) , T = ±, (2.12) 
evM\^Z,a^l'^)s L = ^' y) ' i = 1, 2, 3 , (2.13) 



IRT ^™ 

where: the averages in the l.h.s are defined as in Eq.(1.16) and following lines; — 
X](T ^x,CT^x,o- ^'^^ Px — 5^0- ^x+(Ji o-^x+5i cr- Moreover, if'x, y are two distinct space-time points, 

e'ipl; , = lim F^;f '^(x; y) , r, r' = ± , (2.14) 

e'voip^^; J^),^ = lim Kff'\,c;y) , r = ± , J = 1,2,3 , (2.15) 



e'vl{Ji;4)sL = J™ K^f'^^i^^y), J,/ = 1,2,3, (2.16) 

The proof of this statement is completely analogous to the one in Appendix B of [11] and will 
not be repeated here. Note that, once that the various correlation functions in Proposition 
2 are known, we can reconstruct the functions K^i^{p) and G'2^'^J^(k, p) in Eqs.(1.19)-(1.21) 
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simply by Fourier transformation (provided that the correlation functions in Proposition 2 
have good enough decay properties so that the Fourier transform is well defined), e.g., 

GC'i;o(k,p) = ^ / rfx/ dy [ c^ze^>^(--y)+^P(— )((p^ + e-^^^Vf+.J;V'x-.^^^^^^ 

P-^ Jil3,L) J{I3,L) J(p,L) 

GC4(k,p)-|7^ / d^f dy f d^e^^^-^^^'^^-^^j^m^r^^^^^^^ 

J{0,L) J{p,L) J{p,L) '^''^ 

similar formulas are valid for K^i^{p). 

B. Renormalization Group 

The naive perturbative expansion in U of Wm,^,l(^, (f>) is affected by infrared divergences 
due to the singularity at the Fermi points of the free propagator S'o(k), see Eq.(1.36). The 
case with A — has been studied in detail in [11], where it has been shown that the 
apparent divergences affecting the naive perturbation series can be cured by proper iterative 
resummations; these allowed us to recast the original expansion into a new convergent 
expansion (uniformly in M, (3, L) involving a sequence of effective parameters Zh,Vh, playing 
the role of effective wave function renormalization and Fermi velocity at momentum scale 
2^, h < 0, relative to the Fermi points. The iterative resummations can be implemented 
by using constructive fermionic Renormalization Group (RG) techniques; a key point of the 
analysis, which make the construction of the ground state possible, is the fact that density- 
density interactions are irrelevant in a RG sense. In this section we review the iterative 
integration scheme used to compute Wm,/3,l(^, 0), with particular emphasis to the novelties 
due to presence of the external field A. 

1. The ultraviolet integration. 

Proceeding as in Sections 3.2 and 3.4 of [11], we decompose the propagators ^(k) into 
sums of two propagators supported in the regions of ko "large" and "small" . The regions 
of ko large and small are defined in terms of the smooth support function x(i); the con- 
stant oo entering in its definition is chosen so that the supports of Xo (^'\/ + 1^ ~ PfIJ ^^"^ 
Xo{^\l ^0 + \^ ~ Pf\^ disjoint (e.g., Oq = 1/3 is fine). We rewrite gi(k) as 

^(k) = ^(--)(k) + ^(^-)(k) (2.18) 
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where, setting p'^ = (0,P^) with uj — ±: 

^("-)(k) = ^(k) - ^(-•)(k) , y(--)(k) = J2 Xo(|k - p^l)^(k) , (2.19) 

UJ = ± 

Defining V{A, ^, 0) = V(*) + (*, 0) + {A, J), we can rewrite 

eWM,,,^{A,^) ^ J p^^^y{A,^,<i>) ^ J p(^^{-.)) J p(^^(«...))e^{^.*(* -)+*'"" ^</') , (2.20) 

where P(d^'(' ''-^) and P(d^'^" "-^) are the Gaussian integrations associated to the propagators 
^('•'■•)(k) and ^("•"•^(k), respectively. Using Eq.(2.20) we can further rewrite the generating 
functional as: 



= / P{d¥'-^-^) exp { J2 -A.(V{A, + n)} = (2.21) 

n>l 



=: e 



n>l 



where: S^^ is the truncated expectation with respect to the propagator ^("•''•^(k); Fo,m 
is a constant; Vm is the effective potential on scale 0; 5m(^, \ 0) collects the terms 
depending on A, generated by the ultraviolet integration; ) is defined in the same way 
as J (sec Eq.(2.9)) with ^ replaced by ^(*-''-). As proved in [11] (see Eq.(3.36) and Lemma 
2 of [11]), the effective potential Vm can be written as 



oo 



vm(*(.-))=$:(^l^,-- e z e [n*£^^.,»,-.*£;:i;», 

n=l <Tl,...,(T„=t4- Plr-->/02Ti=l,2 kl,...,k2n j=l 

n 

•W^M;2n;p(kl, . . . , k2„_i) 5(^(k2,_i - ks,)) , (2.22) 

where p = (pi, . . . , p2n) and we used the notation 

6{k) = 6{k)6{ko) , 6{k)=L' h,n,G,^n.G,^ W = /^^o.O , (2-23) 

with Gi,G'2 the basis of A* defined after Eq.(1.3). Moreover (see Lemma 2 of [11]), the 
constant Fo,m and the kernels WM;2n;p are given by power series in U, convergent in the 
complex disc \U\ < Uq, for Uq small enough and independent of f3,L,M; after Fourier 
transform, the x-space counterparts of the kernels WM;2n;p satisfy the following bounds: 



/ 



dXi---cix2„[ n l|Xi-X,r-]|W^M;2n;p(xi,...,X2„)| </3|A|(:7-|C/r-^^'"-l> , 

(2.24) 



l<i<j<2n 
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for some constant > 0, where m — X]i<j<j<2n '^^^ limits Fq — limjvf_^oo -^o,m and 
1^2n;p(xi, . . . , X2n) = limM-700 l^M;2n;p(xi, . . . , X2„) exist and are reached uniformly in M, so 
that, in particular, the limiting functions are analytic in the same domain \U\ < Uq. 
The functional Bm{A, v]>(*-^-)^ admits very similar representation and bounds, i.e.. 



m 



n,m=0 ^ £,P k,p j=l i=l 

hi 

n m 

-Y.P^), (2.25) 

j=l i=l 

where a := (ai, . . . , o-„), ^ := (jji, . . . , tt„J (with jli G {+, -, 1, 2, 3}), 7 = (71, . . . , 72n) (with 
7j G {ea;t, int}), k = (ki, . . . , k2n), p = (pi, . . . , Pm) and the field (f can be either or 0, 

depending on the label 7, i.e., ^t.a,ext = <^k,a and 0t,a,int = ^k,^^^- The kernels W^M;2n,m;p,|;7 
are analytic in [/, they admit bounds analogous to Eq.(2.24), uniformly in 7 and M, and 
converge uniformly as M — >■ 00 to limiting kernels denoted by W2n,m;p,ti;7 (^°^ ^^^^ reason, 
the label M will not play any important role in the following and will be dropped from now 
on). The proof of these claim goes along the same hues as the proof of Lemma 2 in [11] and 
will not repeated here. 

The kernels W2n,m;p,i;^ satisfy a number of symmetry properties inherited from the 
symmetries of the action and of the Gaussian integration, described and proved in Ap- 
pendix B. In particular, thinking of W2n,m;i;'y ^ tensors with entries Vt^2n,m;p,|i;7 defin- 

ing W2A^') - W^2,0;(mt,mt)(k' + p<^,k' + p'^), - W^0,2;«,b(p, -p), W^B,a,(k', p) : = 

W^2,l;«;{mt,mt)(k' + p^^ + p, k' + p%, p), 

W20^') = -izoK + ^iikW2 + ujk'2ai) + 0(|kf ) , (2.26) 

Wr,r'{p) = e-'l^'^^-^'\a + aVr') + Odp^ , (2.27) 
Wrjip) = bTpo + 0{\pf), (2.28) 

WjAp) = e'^"(^^-V)(c5,,,, + c') + O(IpI^) , (2.29) 
WrA^', p) = Ao + \3ras + 0(|k'| + IpI) , (2.30) 
Wj,Uk',p) = Aie-^(^-^)-3^^ + o(|k'| + IpI) , (2.31) 



where the constants zq, zi,a,a' ,b,c,c' , \q, Xi, X3 are all real. The proof of Eq.(2.26) is in 
Lemma 3 of [11], while the proof of the remaining relations is in Appendix C. 
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2. The infrared integration. 



After the integration of the ultraviolet modes, we decompose the infrared propagator as 
a sum of two quasi-particle propagators: 



a)=± 



(2.32) 



— * 

where, if k' = (/cq, k'), 



k'^, \-Von{k'+p^) -iko 



(2.33) 



andBl^ ^ {2n/3-\Z+^)n{ko : Xo{2-^\ko\) > 0}} x{fGi + fG2-p^ , 0<ni,n2<L}. 
Correspondingly, we rewrite ) and 0^ as sums of two independent Grassmann fields, 
"^ia^^ = J2u:=± e'^'^'^'^^w'^ and = Z)a;=± ^'^''''^''^x.a.a;' ^^d rewrite the generating 
functional as (dropping systematically the M label): 



(2.34) 



g-/3L2i.o+5(>o)(^^^) J p^^^^^^^(<o)^gV(o)(*(^o))+B(o)(A,*(<o)>)+(A*(<o),Ti*(<o))+(<;!,,*(<o)) ^ 



where: 



Xo(|k' 



)>0 



11 11 "^k',a,p,a;"^k',a,p,a, 



(2.35) 



a=nX0{\^'\)>0 

■ exp { - (/3L^)- xo-^(ik'i)^g°!:a,.(ko^gt ) 



and 



Co,a.(k') 



-iko -voVL*{k' +p^\ 
-Von{k'+pp) -iko 

—iko vo{ik[ — ujk'2) 
Vo{—ik'i—(jjk'^ —iko 



(2.36) 



;i + i?o,.(k')). 



with J\fo chosen in such a way that J P;^g,Co(c?^*^~°^) — 1, Zq — 1 and Ro,oj a matrix such 
that ||i?o,u;(k')|| < C|k'|2. Moreover, the functionals V^^\ (A*(^°), Ti^^^^)) and (0,*^-°)) 
are the same as Vm, (A, J^^ ""-)) and (0, ^f^^-^)) in Eq.(2.21), with and 0± rewritten 
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in terms of ^^L""^ and 0^; similarly S^-^\A,(f)) + 5^°^ is the same as Bm after the same 
rewriting, with 5"*^-°^ collecting all the terms depending on A and/or (f) but not on ^(^°). 
E.g., to be more explicit, defining ^^7^^ := A^u: _^u:' a and recalling the definition of n±, 

P )R Pf Pf "tP li* 

i.e., n± = (1 ±(j3)/2, 

(^^(<o), = -i- ^ 5] i^rr'^gi;,,.r-1k', pO^gt" , (2.37) 

^ u;,u>',a,t|k',p' 

with r;;f tk', pO = rir and 

Tj^f (k', p') =vo- e^¥0--i)(-'-+-"-) ^^^^-ik'(5,-5i) _ ^_g+i(k'+p')(5,-50) . (2.38) 

' 3 

similar expressions are valid for the other functional appearing in the exponent of Eq.(2.34). 



The integration of Eq.(2.34) is performed in an iterative fashion, justified by the in- 
frared divergence caused by the singularity at k' = of the propagator ^ir°^(k'), see 
Eqs.(2.33),(2.36). We define //,(k') := xo(2"^1k'|) - xo(2"''+^|k'|) and rewrite 

Xo(lk'l) = /'^(kO X.(k') + A+i(k') + • • • + /o(k') . (2.39) 

fc=— oo 

The integration procedure consists in the following: we first rewrite Xo = X-i + /o and corre- 
spondingly decompose the propagator g^^r"^ = 5'ir~^^+5'(?^ and the field = ^ir~^^-|-^L°^; 
next, we integrate out the "0-mode" , i.e., the Grassmann field and re-express the result 
in terms of a new effective potential on scale —1; then we iterate, by integrating out step by 
step the degrees of freedom on scale —1, —2, . . . ,h + 1, with h < 0. After the integration of 
the fields on scales —1, —2, . . . ,h+l, we get: 

where Fh, Ch, V^'^\ B^'^\ Th, Q^^'> will be defined recursively (in particular, Fh, Ch 
[called Ah in [11]], V^'^), Q^'^) have already been defined in [11]) and P^^^cM-^^-^^) is 
defined in the same way as P^^^Coid^'^-^'') with xo, C'o,^, -^o, 'J^o, -Ro,w replaced by 

Xft, Ch,uj, Zh, Vh, Rh,uj, respectively Moreover, 



oo n 

n=l o;,P,'id,k' J=l 

2n 

• ■ • > k^n-i) 5(E(-1)'(P^^ + K)) (2-41) 
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and 



V ' ' ^/ / ^ [ 11 ^k2^_l,0-3,p2j_l,72j-l,W2j-l^k2^-,0-3,P23,72j,W23 

n,m>0 <I>P>tt>7)'t^ i=l 

2n 



(2.42) 



1=1 j=l i=l 



int 



where: the *'s on the sums remind the fact that only terms exphcitly depending on (p^, 
^(</i)± contribute to 5*^'*^; ui — {ui, . . . , U2n] oJi — oj[, ■ ■ ■ , cOm — co!^)- 
The terms (^*(^'^), and (0, Q('^+i)*(^'^)) read: 

k',p' 



w.CT.k' 



for suitable matrices Tj"^'^^(k', p') and Qk't)^''^' "^^^^ ^k''*w ^ iterative procedure 

goes on up to scale /i^j, where /i^j is the largest scale such that 007'*''"^ < ^. The result of 
the last iteration is e^^^^^^'^). 



3. Localization and renormalization. 

In order to inductively prove Eq.(2.40), we rewrite 

yW(^{<h)-) ^ £,\;{h)(^q,{<h)^ ^ Ji\;ih)(^q,i<h)^ ^ (2.45) 

where 

£V('-)(fl>«'") = ^E E*S!>5U,(k')*gt . (2.46) 

and TZV'^^^ is given by (2.41) with J^^i replaced by X]^2' ^^^^ contains only the 
monomials with four or more fields (note that in Eq.(2.46) the two fermionic fields have the 
same a;-index; terms with two different quasi-particles indices are not allowed by momentum 
conservation, see also the remark after (3.62) of [11]). Moreover, defining 

<!.'(k', pO := E ^S«;(.,.';.-.');M,n*)(k' + p', k', p') , (2.47) 

Wl,W2 
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we rewrite: 

B'^^\A,^^^^\(P) = CB'^^\A,¥^''\(P) + nB^''\A,¥^''\(P) , (2.48) 

where 

+ E E p')*gt!: . (2.49) 

^ a;,a;' k',p' 

At this point we "reabsorbe" CV^^^ in the fermionic gaussian integration and CB^^^ into the 
definition of the effective source terms: 

^ g-/3L2(Fh+eft)+5(^'')(A,<^) j ^ ((^^(<M)e7^vW(4'(^'')) . ^2.50) 

7es('')(yl,*(^''),</-)+(A<I'(^''),Th1'(=^''))+(0,Q('')*(^'')) 

where Ch is a suitable constant (see Eq.(3.67) of [11]) and 

C,_i,^(kO = C,,^(k') +X/.(k')W^2;(L,.)(k') , 

fc=/i+i 

«K = qK"* + ( E Qi.'^si")*'^-)*'''' • 

fc=/l+l 

T-^-'(k', pO = r„;^^^^;(k', pO + w„fj,^,(k', p') . 

The second and third equations in (2.51) can be proved as in [11], see Eqs.(3.111)-(3.113). 
We are now ready to perform the integration of the ^^'^-^ field: we rewrite the Grassmann 
field as a sum of two independent Grassmann fields + and correspondingly 

we rewrite the r.h.s. of Eq.(2.50) as 

^-0LHF,+e,)+si>^HA,^)J p^^_^^^^_^(d*(<'^-i))yp^^^^^_^(d*W)e^v('')(*(^''-)+*('^)) . 

.g7eB('')(A,*(^''-i)+*(''),0)+(A(«r(<''-i)+*(h)),Th(*(^''-i)+*('')))+(0,QW(*(^''-i)+*(''))) 52) 

where 

a-i,u,(k') = a,.(k') + w^2;(i,,)(kO . (2.53) 
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In [11] we proved (see Eqs.(3.68)-(3.69) and proof of Theorem 2 of [11]) that the single scale 
propagator, defined as 

£'^)(xi-X2) = -^ J2 e-'''(— ^)/,(kO[Z,_i,.(k')]"\ (2.54) 
can be rewritten as 

^W(xi-X2) = -^ J2 e-^'(---)- (2.55) 



-1 



1 + Rh-i,o.{k'))-' , 



Zh-i(k') \vh-i(li'){-ik[ - uk'^) -iko 
with Zh-i,Vh^i two functions such that (choosing < ^ < 1) 

|Zft_i(k') - Zh{k')\ < (const.)|t/|2^^ , Zo(k') = Zo = 1 , 

\vh-i{k') - %(k')|| < (const.)|C/|2^'^ , vo{k') ^Vo^^t (2.56) 

and Rh-i,uj a matrix such that \\Rh-i,oj{k')\\ < (const.) |k'p and, if 2^'^ < |k'| < 2^*+^ 
||9k'-Rfe-i,u;(k')|| < (const.)2(2-")'^. The constants := ^/^(O) and t;^ := %(0) play the 
role of effective wave function renormalization and Fermi velocity on scale h. Similarly, the 
local parts of the matrices T^^j^ (k', p') play the role of effective vertices; in particular, the 
"relativistic" vertex functions, which represent the dominant contribution in the infrared 
to the kernel of the three point function, are defined as (see the subsection The three-point 
function below for more details) : 

3 3 

Zo,H = ^rT(0' 0) ' = E(^^)i^r(0, 0) , Z,,, = -ua, $^(^,)27Xr (0, 0) . 

T=± j = l j = l 

(2.57) 

From the symmetries discussed in Appendix B and C (see also Eqs.(2.30)-(2.31)), Z^^, 
H = 0,1,2, are all real, independent of u and proportional to the identity matrix and, 
therefore, they can be regarded as constants, as we will do in the following; moreover, 
Zi,h = Z2,h, Zo,o = 1, ^1,0 = Z2,o = vo and \Z^,h-i - Z^,h\ < (const.) |C/|2^^ see Eq.(2.64) 
and following discussion for a proof. 
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Now, going back to Eq.(2.52), if we integrate out the field on scale h and define: 

^ g5(S:'')(A,</>) J Pf ^^^(?i)-jg7ev('')(*(^''-i)+*(''))+7^BW(A,*(^''-i)+*(''),<^) _ 

(^^(h),T;^$(<'>-i))+(A*(<''-i),r,,*(''))+(A*('') ,rft*(''))+(fli,Q('')*('')) 

wc get Eq.(2.40) with h replaced hj h — 1 (and Fh-i = + + eh). 

The integration in Eq.(2.58) can be performed by expanding in series the exponen- 
tial in the r.h.s. and integrating term by term with respect to the gaussian integration 
Pf^^^_^{d^^^^). This gives rise to an expansion for W/3,l(>1, 0), which can be conveniently 
represented in terms of Gallavotti-Nicold trees, as described in Section 3.3 of [11] and in the 
next subsection. 



4- Tree expansion. 

For each n > and m > 2, we introduce a family 7^"^ of rooted labelled trees, defined in 
a way similar to the family Th,n described in Section 3.3. of [11] (which we refer to for more 
details), with the following minor modifications: 

1. 7^™ has n+m endpoints (rather than n); n of them are called normal endpoints and m 
of them are called special endpoints; moreover, the special endpoints can be either of 
type A or of type 0. If is a normal endpoint, then, as in [11], it is associated to one of 
the monomials with four or more Grassmann fields contributing to 7?.V^°^(^'^-'*''~^^); 
if V is a special endpoint of type A, then it is associated to one of the monomials 
contributing to (A^^^^^--^), T;,„_i^(^''''-i)) - (A^^^'^''-^), T/,^_i^(^'^''-2)); if t; is a spe- 
cial endpoint of type 0, then it is associated to one of the monomials contributing to 
(0,Q('^''-i)*('^''-^)). 

2. Among the sets of field labels and external field labels associated to the vertex v, 
denoted by and P„ in [11], we distinguish the field labels of type A, ip and 0; 
we correspondingly introduce the sets P/, P^, P^, If and I^. All the trees 
contributing to Wj3,L{A,(f)) + ^L'^Fhp are characterized by the fact that P^ = and 

Vn ^ vn r V. 
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Given r e 7^^, the sets Pf, # e {^4, 0, ■0}, satisfy several constraints that depend on r 
(see [11]). In particular: 

a. denoting by Vq the vertex immediately following the root on r, we have |P^| = and 
\Pvo\ + \Pt\ > 0; if ^; > v^, then \Pt\ > 0; 

b. for any e r, sets of external field labels such that \Pv\ = \P^\ — 2 are not allowed 
(this follows from the definition of localization Eqs.(2.45)-(2.46) and from the choice of 
"reabsorbing" at each step the local part into the gaussian integration, see Eqs.(2.50), 
(2.52), (2.58)); 

c. if V is not an endpoint, then sets of external field labels such that |P„| = 2 and 
\P^\ = \p^\ = 1 or such that \P^\ = 3, |P/| = 1 and \P^\ = 2 are not allowed 
(this follows from the definition of localization Eqs.(2.48)-(2.49) and from the choice 

of "reabsorbing" at each step the bilinear terms (0, W^2l'^) i^^*^ definition of 
(0,Q(^)^) and the "vertex" terms (A^, VT^'^'^) into the definition of (A^,Tft_i^), 
see Eq.(2.51)). 

As in [11], we denote by Vr the family of all the choices of P* compatible with these 
constraints and by P the elements of Vr- The generating functional can be expressed as a 
sum over trees in the following fashion (analogous to Eqs. (3. 77), (3. 79), (3. 87), (3. 88) of [11]): 

>v^,^(A0)+/3^X = EE E E E E^^'^(^'P'^)' (2-59) 

n>0 m>2 h>h^ t^'T^^„ PsPt TeT 

where, as explained in [11], T is a suitable family of spanning trees. The contribution 
VV^'*^(t, P, T) can be further rewritten as 

W('^)(r, P, 3^) = / ^x„„A(P,^) 0(P,t) <^^,t(x.o) , (2.60) 

where 

^y^vo) - 11 e ^x(/),«(/) ' n^vo) - 11 e y^^{f),a{f)Mf) 

(2.61) 
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and, calling v*,. . .,v* the endpoints of r, putting hi — h^* and denoting by K^'i'^'hcy*) the 
kernels associated to the endpoints, 



W^t,P,t(x„o) = 



,i=l 



(2.62) 



J] i- /dPT.(t.)detG''^'"^''(t,) 



V not e.p 



In the latter equation, dPr(t) is a probability measure with support on a set of t such that 
tjj/ = Uj ■ Uj' for some family of vectors u, e of unit norm. Finally G^''^{t) is a Gram 
matrix, whose elements are given by, see [11], Eq.(3.83): 



Pi 



(2.63) 



this matrix takes into account all the possible contractions of fields not involved in the 
spanning tree T. See [11] for more details. The effective potential V^'*-* and the effective 
source term S*^'*) admit representations very similar to Eqs.(2.59)-(2.60)-(2.62), the main 
difference being that these are expressed as sums over trees and field labels satisfying slightly 
different constraints: i.e., the trees contributing to V^'*^ do not have special endpoints and 
are associated to external field labels such that = P^^^ = and |P^J > 0; the trees 
contributing to B^^^ have at least one special endpoint, and are associated to external field 
labels such that \P^J > and |Pi| + |P„*J > 0. 



5. The kernels of the special endpoints of type A 

In order to prove Theorem 2 and Proposition 1, we will be particularly concerned with 
estimating the contributions with \P^J =2 and \PfJ = or with \P^J = 1 and |P/J + \PfJ = 
2. The key estimate that we preliminarily need to prove is 

ll<j,.'(k^p')ll<(const.)|C/|2^^ (2.64) 

for all /c < and with 9 e (0,1), uniformly in k',p'. Note that Eq.(2.64) implies, in 
particular, that the kernel of the special endpoints of type A is uniformly bounded as 

\\T-/{k',p')\\<Co, (2.65) 

for all /c < and a suitable constant Co, and that iZ^^h-i — Z^^h\ < (const. )]C/]2^'', as claimed 
after Eq.(2.57). Let us proceed by induction: we assume the validity of Eq.(2.64) ior k > h 
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(so that Eq.(2.65) is valid for all k > h) and prove it for k — h. Using the tree expansion, 
we can rewrite: 

P n>l r&T} PePrTeT-^ 

(2.66) 

where the * on the sum reminds that P^,,, = U P^, = {/i} ^-^d P^^^ = {/2, /s} (with 
uj{h) = oj{f2) = CO, cu'ih) = cuifs) = u, x(/i) = z, x(/2) = X, x(/3) = y, tt(/i) = i 
£(/2) = "^(/s) = ~ s-iid '^(/2) — (^(/s))- Using translation invariance, the representation 
Eq.(2.62) and proceeding as in the proof of Theorem 2 of [11], we get (see Eq.(3.93) of [11]) 



0',P')ii<E E EE/ U'li-'-y:'' 

n>i reri FeVr tgt zeT* 



UK- 



(hi) I 



i=l 



(2.67) 



Yl max |det G"^-^" (t.) | H I kS"^ " 



sJ. tv 

_v not e.p. lETy 

The r.h.s. of this equation can be bounded dimensionally, using the scaling properties of 
the propagators ^(^^"■'(x) and of the Gram determinants det G'*'"'^''(t^) (see Eqs.(3.92),(3.94) 
of [11]). Following the proof of Theorem 2 of [11] and using in particular Eqs.(3.94)-(3.95)- 
(3.96), we get the analogue of Eq.(3.97) of [11], that is 

IK^l(k',p')ll< 



(2.68) 



SEE EEc'in^^ 



hv{{j:Zi\Pn\)-\Pv\-3{sv 



V not 
e.p. 



-■''][n 



C 2 



\p,f\ 



V e.p. 



where Pi — \Py* \ and C is a suitable positive constant, larger than the constant Co appearing 
in Eq.(2.65) (in fact, in deriving this bound, we estimated the kernel of the special endpoint 
of type A by using the inductive hypothesis Eq.(2.65)). The r.h.s. of this expression can 
be rewritten in a convenient form, using the analogues of Eqs.(3.98)-(3.100) of [11], that is 
(recalling that m is the number of special endpoints in r - equal to 1 in the current case) , 

E [(E i^^i) - = Mi^Si - inti) + E (^^ - - ml) , 

h not i=l h not 

e.p. e.p. 

E hv{sy - 1) = /i(n + m - 1) + E ~ hy,){n{v) + m{v) - 1) , (2.69) 



h not 
e.p. 



h not 
e.p. 



where: v' is the vertex immediately preceding v on r; 1^ D Py is the set of field labels 
associated to v (i.e., including both the internal and the external fields to v); n{v) is the 
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number of normal endpoints following v on r; m{v) is the number of special endpoints 
following V on r. Note that in the current case, where there is only one special endpoint of 

type A, ni{v) = \P.^\ and m = l-P^I- Plugging Eq.(2.69) into Eq.(2.68) and using Eq.(3.100) 
of [11], we get the analogue of Eq.(3.101) of [11], that is. 



;k',P')ll ^ XI XI ^C''^2''(3-I<I-I-p4I) . J JJ^ l_2^K-\,)(3-\p*\-\p,^\)'^ . 

n>l rer.i PeVrTeT vnot^'"' 

~ ^'h,n e.p. 

.[ Yl 2'^'''(l^'l+l^"l-^)] [ n C^|[/|^-i] , (2.70) 

V e.p. V e.p. 

where we used that, if v is an endpoint, then = \P^\. At this point, by Eqs. (3.102)- 
(3.103) of [11] and by the argument described after Eq.(3.103) of [11], we end up with the 
analogue of Eq.(3.104) of [11], that is 

n>l 

Recalling that in the current case 3 — |P^J — = 0, this proves the desired estimate on 
the kernels of the special endpoints of type A, Eq.(2.64). A similar strategy also allows us 
to prove that 

7^;:;;^'(k',p') = t;;;;-'(o,o) + o((|k'| + |p'|) + \u\ {\k'\ + \p'\f) , (2.72) 

with 9 e (0, 1). Let us also recall that also the kernels Q^^^ of the special endpoints of type 
admit a uniform bound of the form < Co, see Eq.(3.114) of [11]. We are now ready 

to give the proof of Theorem 2 and Proposition 1. 



6. The three-point function (Proof of Theorem 2). 

The goal is to bound G2'^/(k, p) at k 7^ p^, p 7^ 0, with |k — p'^| and |p| sufficiently 
small (and |p| <^ |k — p"^!), for a given u e {+, — }. The three-point function, by definition, 
using Eqs.(2.10)-(2.13)-(2.17) and Eqs.(2.59)-(2.62), can be rewritten as 

G&(k,p) = 3^Efe)4E E E E E Ax„e*(--''H*(-)w.-«,(x„), 

^ j=i n>o h>hp TeT^^PeVrTeT-^ 

(2.73) 

where the ** on the sum over P reminds that P„o = P^UP^^, — {fi} and = {/2, /s}, 
with uj{fi) = uj'{fi) = u{f2) = ^(/s) = t^, x(/i) = z, x(/2) = X, x(/3) = y, iJ(/i) = j. 
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£(/2) = "^(/a) — ~ and cr{f2) — cr{f3). The reason why all the quasi-particle indices of 
the external legs are equal to u is that, by assumption, |p| <^ |k — p"^] <^ 1, so that by 
momentum conservation all other choices of quasi-particle indices give zero contribution to 
G'2']^j(k, p). The trees contributing to G'2']^^(k, p) have a few more features that are worth 
remarking. First of all, among the three special endpoints of r £ 7^^„, one of them is of 
type A (let us call it va and note that /i e P^) and the other two are of type 4> (let us call 
them vfjVT, with /2 G P*^. and fs e P%)- Moreover, let be the (negative) integer such 
that 2^^ < |k - p"^! < and let |p| be so small that 2'^^-! < |k + p _ p<^| < 2^^+'^: 

then, only trees with \h^± — /ik| < 1 and with /i < /ik + 1 contribute to G2'H(k, p). 

Now, let us distinguish in the r.h.s. of Eq.(2.73), the contributions with n = and those 
with n>l. The latter can be bounded in a way completely analogous to w!^!"^ ^ p'), with 
the following important differences. In the current case m = 3 and m{v) — m^{v) + m'^(i'), 
with ■m'^{v) — \Py\ (resp. m'^{v) — \P^\) the number of special endpoints of type A (resp. 
type (p) following v on t. Taking this into account and following the same strategy used to 
bound W^^^cj^ Eq.(2.67)-(2.68)-(2.69), we get the analogue of Eq.(2.70), that is we can 
bound the contributions with n > 1 to the three point function by 

n>l h<hi^+l tGT? PeVrTeT vnot^'"' 

e.p. 

.J JJ 2'*"'(l-P'^l+l-f'^l+l^'l-3)] [ Y[ , (2.74) 

V e.p. V e.p. 

Now note that, if is a special endpoint of type 0, then \P^\ + \P^\ + \P^\ — 3 = — 1 and 

that, if ^; is a special endpoint of type A, then \P^\ + \P^\ + |P^| —3 = 0. Therefore, using 

the fact that \Pvo\ — 3, we can rewrite Eq.(2.74) as 

** 

g2-2'»k^ XI XI xi^"' [ n — j2('^''-'*''')(^-i-^''i)] • 

n>l h<hk+l reT? PeVr TeT v not 

- - '^'h,n e.p. 

.[ 2''^'(l^''l-3)][ J] C^lt/l^-i] . (2.75) 



V normal v e.p. 

e.p. 



The only potentially dangerous contributions to Eq.(2.75) are those coming from a vertex v 
that is not an endpoint and such that |P^| = 3. By construction, such vertex has necessarily 
|P/| = \Pf\ — \P^\ — 1; on the other hand, by momentum conservation, < /i„ — < 2, 
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simply because |p| <^ |k — p"^! <^ 1 and, therefore, the quasi-momenta associated to the 
external fields and ^ have essentially the same momentum scale (i.e., the two scales differ 
at most by 1); in conclusion, the overall contribution coming from a vertex that is not an 
endpoint and such that \P^\ = 3 can be bounded by an 0(1) constant and gives no trouble. 
Therefore, Eq.(2.75) implies that the overall contribution to G2 H(k, p) coming from trees 
with n > 1 can be bounded by (const.) |t/|2(-2+^)''''. We are left with the contributions 
coming from the trees with n — 0, which read (defining k' = k — p'^) 

e E Qa;,.^i'^(k'+P)[E(5■)^^?^;:r(k^p)]£'^'^k^^ (2.76) 

h,h'=h]^ j=l 

where h :— max{/i, /i'}. Using Eq.(2.64), as well as the estimates on Q^'^^ and on the 
two-point Schwinger function proved in [11], see Eqs.(3.114),(3.120)-(3.122) of [11], we can 
rewrite Eq.(2.76) as 

e (^(k + P) [E(^^)^^^rjk', P)] S{k)) (l + 0{\U\2''^-)) , (2.77) 

for some 9 e (0, 1). Moreover, using Eq.(2.72), we can further rewrite Eq.(2.77) as 

e (^(k + p)Zi,hJiip^, 0)5(k)) (l + 0(2'^^) + 0(1^712'^'''')) , (2.78) 

where the vertex functions ^ were defined in Eq.(2.57) and we remind the reader that 
Ti{pp, 0) = — (J2 and T2{pp, 0) = —coai. Theorem 2 is an immediate corollary of the previous 
estimates and, in particular, of Eqs.(2.75)-(2.78). Using the fact that \Zi^h — Zi^h-i\ < 
(const.) 1 1/ 1 2^^, we also find that the constants in the statement of Theorem 2 coincide 
with the infrared limit of the running coupling constants Z^^h, i-e., Z^ — lim/i_^_oo Z^^^h- ■ 



7. The response function (Proof of Proposition 1). 

In order to prove Proposition 1, we start by deriving bounds on the current-current 

response function Kj^^^'x — y) at distinct space-time points, x 7^ y, which can be expressed 
in terms of the tree expansion as follows: 




(2.79) 
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where the * * * on the sum over P reminds that P^^ — — {/i, /2}, with uj{fi) — uj'{fi) — 
^{f2) - uj'{f2) = 0, x(/i) = X, x(/2) = y, x(/3) = y, = j and ^(/a) = /; moreover, the 
* over c?x^Q reminds that we are integrating over all the variables in x^^ but x and y, which 
arc fixed (and distinct); finally, Hjj/ is the Fourier transform of Titrj-i^ and {Hjj/ * W^py) 
denotes the convolutions between Hjj' and H^t,p,t- The trees contributing to K^^^^^ have two 
special endpoints of type A; we call them vi and V2, we denote by v* the rightest vertex such 
that V* < Vi,V2 with respect to the partial order induced by the tree and by h* its scale. 
Proceeding as in the previous subsection, we distinguish the contributions to coming 
from trees of order n = or n > 1; we denote the two by K^^ and K^^, respectively. The 
latter can be bounded by (using a notation analogous to the one used in Eq.(2.67)) 

y)ii < E ii^-'ii^. E E E E E / ri di^i - ■ (^-so) 

j,j'=l n>l h>hi3 T^Tj^^PeVrT^T'^ leT* 



i=l 



n ^max|detG'^'"^^(t,)| J] M'^^i - Yi] 

. V not e.p. lETy 

where Yii&T* product over all fines of the (modified) spanning tree T*, but one fine 

belonging to the subtree Tx^y C T* connecting x with y and contained into the cluster v* 
but not in any smaller one; let us call I this special line. Eq.(2.80) can be bounded in a way 
analogous to Eq.(2.67), with the important difference that the Li norm of the propagator 
associated to I is replaced by its norm, which has a factor 2^'^* more as compared to 
the Li norm; moreover, in order to take into account the decay between x and y, we can 



extract a factor 



Cn 
A 



i+(2'i*|x-y|p from the product of the propagators in the spanning tree (here 
N > 1 and Cjv is a suitable positive constant); we are still left with an expression that can 
be estimated in the same way as Eq.(2.67), thus leading to the upper bound 



l<Hx-y) 



'oEE E EE^"2M3-ip.oi) 

>ih>hp reT^^P^VrTeT 



2ihy-K'){3-\Pv\) 



V not 
e.p. 



c^\u\ 



-,3h* 



l + (2'»>-y|)^ 

(2.81) 



V e.p. 



V e.p. 



Proceeding as in the previous subsections, we see that this can be further bounded as 

, , nn 



|i^^^(x-y)|| < ^ ^ ClC/r2'^(^+^)2 

n>l h>h(i h*>h 



l + (2^*|x-y| 



liV 



(2.82) 
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Now, noting that, by construction, h > h* and exchanging the order of summation over h 
and h*, we get (picking N — 5 and \U\ small enough) 

ll<Hx-y)|| < (const.)e^|C/| Yl i + (2^'U-^l)^ ' (^.83) 

h*>hp ^ ' 



which implies 

ll<Hx-y)|| < iconst.y\U\^^j-L-— (2.84) 

and, therefore, in the limit P,L ^ oo (which exists by the uniformity of the bounds and by 
the term by term convergence of the series, see Appendix D of [11] for more details), the 
Fourier transform of K^^^ is continuous and continuously differentiable for all p e M x B (in 
particular at p = 0) . 

We are now left with the contributions to the response function coming from the trees 
with n = 0, which read 

3 



j,j'=l h,h'>hB uj,Lo' ^ ' ' ' 



■"k-i^.f "-«-*>-Ap-)'^'«""(''' + p'^^J'*"' p')£';''(k')r;,^(k' + p'. -p')} . 

where h — max{/i, h'} and the factor 2 in the r.h.s. takes into account the summation 
over the spin degrees of freedom. It is important to notice that the integral in the latter 

expression can be rewritten as c'^^^f-Pf ■"•^~^''/)^'^ (x — y), with (x — y) a function 

admitting the dimensional bound: 

nn-max{h,h'} 

|5x/r,.^'-'(x)| < ^^2^^^"'' \ + (2n^ax(M-}|x|)^^ ' ^'-''^ 

for all N > 0. Eq.(2.85) implies that the contributions to K^^\-x. — y) with uj ^ uJ can be 
rewritten as e^^'^'F~Pp'^^^~'^^F^^_^^{yi — y), with 

onmax{/i,/i'} „/ 

|a^F,_,(x)| < e^c^y 22(''+''') = jTTJxT^ < tht— , (2-86) 

h,h' ^ I I'' II 

for all A?" > 4 + n. This implies, in particular, that we can rewrite 

3 



_pip(^-y)rr' rr' 
2 ^ 'Ip 'l-f 



if<°'(x-y) = E (m-U EE/ 

j,j'=l h,h'>h0 OJ ^ ^1 

■Tr{9!^\^' + p)T;i'i^', P)e'nk')7^;'f (k' + P, -P)} + Hl^i^ - y) , 
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where the Fourier transform of H^^\x — y) is continuously differentiable in a neighborhood 
of p = 0. Using Eq.(2.72) and the definition of Zi^h, we can further rewrite this expression 
as 

h,h'>hp uj J \ ) \ \ 

•TV{e)(k' + p)r,(pT, o)£'^')(k')r^(pT, 0)} + :ffS(x - y) , 

where the Fourier transform of i?*^"^ (x — y) is continuously differentiable in a neighborhood 
of p = 0. Finally, rewriting Zi^^ — Zi -\- 0{\U\2^^), using the expression of the two-point 
function in terms of a sum of single scale propagators (see Eq.(3.121) of [11]) and using the 
definition of Vik^p) in Eq.(1.39), we get 

+^^Hx-y), (2.87) 

where the Fourier transform of H^^\yi — y) is continuously differentiable in a neighborhood 
of p = 0. Combining Eq.(2.87) with the explicit expression Eq.(1.42) for the two point 
function and with the bounds derived above on K^^\ we finally obtain the statement of 
Proposition 1 (the parity property of -R(p) stated in item 3 of Proposition 1 easily follows 
from the symmetry properties listed in Appendix B and C). ■ 

Appendix A: The conductivity of the non interacting system 

In this Appendix we prove Eq.(1.40), with aij\u=Q given by Eq.(1.38): 



^3 n po^o+ J In Jq \b\ Po 

(A.l) 



y|c7=o 3^ ^ 27r 
where \B\ = 87rV(3y3) and 



, ai(k)\ 

m,0)= . M , (A.2) 
am ) 

with ai{k) = f [1 - e^i*^i cos {^k2)] and a2{k) = ^e^i*^i sin {^k2). 

Let £ > be a small number, independent of po, to be eventually sent to zero. In the 

— * 

integral to be evaluated, we distinguish between the region where \^{k)\ > e and the region 
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where \^{k)\ < e: 

■So{\^)V,{k,Q)]-[xmk)\>e) + xm^)\<^)\ , (A.3) 

The integral associated to the region x(|Q(^)| > e) is uniformly convergent as po — >■ 0+: 
therefore, we can exchange the integral with the limit and check that the integral of the 
limit is zero (simply because the integrand is odd in kg). Next, in the integral associated 
to the region < e) wc rewrite the propagator as the relativistic propagator plus a 

correction (similarly, we rewrite Fj as its relativistic limit plus a correction). The corrections 
are associated to absolutely convergent integrals, uniformly in po as po ^ 0"*", and one can 
easily check that their contribution after having taken £ — > is equal to zero. We are left 
with (after having changed variables and having included a further factor 2 coming from the 
summation over the two Fermi points): 



TT /i £^Opo^O+ I^PO J '^7^ J\k'\<e k^ + \k'\^ ^ {-ik'^-k'^ ik^ 



1 (iiko+po) ik[-k'2 \ 1 I iko ik'^-k'^ 

{ko+PoY + \-ik[ - k'^ i{ko+po) j k^ + \k'\^ y^tk[ - ^ tko 
where (Tj, i — 1,2, are the first two Pauli matrices. Now, if i 7^ j, the r.h.s. of this equation 
is equal to 

±-— hm hm — / — ^ / dk! ^ —] , (A.4) 

n h e^Op,^o+poJ 27r J\k'\<e k^ + \k'\^^{ko+poy + \k'\^ kl+\k'\^^ 

— * — * 

which is zero by the symmetry under the exchange k'i — > — k'. If, on the contrary, i = j, 
we get: 

I = lim lim 1 / ^ / dk'^— \( - _JMh±P2l_] + 

\U=0 TT h e^Opo^Q+poJ 27r 7|fe,|<, k^ + \k'\^ ^k^+\k'\^ ( kn + Pn)^ + \k' \^ 



ki + \kf {ko+Po)^ + \k' 

(A.5) 



+{-iy[2{k[y - 2{k'^y] ( — U — ^ 

^ki + \kf {ko+Poy + \kT' 
Now, the terms in the integral proportional to 2{k[y — 2{k'^'^ are zero by the symmetry 
under the exchange k'-^i — )-/c2- Therefore, we are left with: 



U-^ TT h e^Opo^O+ Po J 271 J\%,\<, k^ + \k'\'^ + \k'\^ {ko + pof + \k' 



(A.6) 
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that is 



(^n rr-n — 1™ 1™ / dk ■ k 

It^-U /l £^Opo-^0+ Po Jo 



f dkp / 
J 27r\ 



ko{ko+po] 



(A.7) 



The integral in ko can be evaluated by residues to give: 



dko 
2wi 



A-2 



ko{ko+po) 



{kl + k^y [{ko+Po? + k^] ■ + 



ftp 



{ko + i/c)^ 



ko=ik 



ko{ko+Po) 



[{ko+PoY + k^] ■ [ko + tk] 



ko=ik 



koiko+po] 



[ko +Po + ik] ■ [k^ + k^] 



ko+po=ik 



that is 



/ 



dko 
27 



^2 



ko{ko+Po) 



Pi 



.{kl + F)2 [{ko + Po)' + A;2] • [kl + J Ak{pl + 4F) " 
Plugging (A. 9) into (A.7) gives 

g2 /"£ 1 e' e' vr 

Cii rr n = 8— lim lim — / dk— 75-— = — lim lim arctan(2e/po) = - — , 



(A.8) 



(A.9) 



(A.IO) 



which is the desired result. 



Appendix B: Symmetry transformations 

In the Appendix we collect some symmetry properties of the fermionic action, i.e., some 
transformation of the fermionic fields and of the external sources that leave separately in- 
variant both the gaussian fermionic integration P(d^) and the interactions V(^), (^,0), 
{A, J). These symmetries will be also preserved by the multiscale integration and, therefore, 
they will allow us to exclude the presence of possibly dangerous terms in the effective action 
at scale /i, see Section 2 and Appendix C. In the following, we denote by 0-1,0-2, 0-3 the stan- 
dard Pauli matrices and we use the foUwoing convention for the Fourier transform of the A 
field: = (/3L2)-i ^p^g^^ e-'PMp,„, where B^^l = 27r/3-iZ x Bl and ^ e {+, -, 1, 2, 3}. 

Lemma B.l For any choice of M,I3,L, the fermionic Gaussian integration P{d'^), the in- 
teraction V(\l/) and the source terms (^, (p), {A, J), defined in Eqs.(2.3)-(2.7), are separately 
invariant under the following transformations (here 7 e {ext, int} and (fi]i,a,int — ^k,a) while 
'f]s.,a,ext — (t>\i,a; whenever this will not create ambiguities, we shall drop the labels a and 7, 
i.e., we shall use (p^ as a shorthand for (p^fj^): 
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(1) Spin flip : (p^^,^ O 'fl,-a,'y; 

(2) Global U{1): (fil,^^^^ ^'^°"'^k,a,'y> ^^^^ eR independent ofk; 



(3) Spin SO{2) : j '^'"'^''^ ) ^ e'''^"^ [ '^'''^''^ VwitheeT^ M/27rZ independent ofk; 

(4) Discrete rotations : (^^ ^ e-*^(^"3-*'"i)^(^^k; ^ (^^^6'^^*"'"^'^^; ^p,± ^ 
^rp,±e^^2(*'"3-«i) ^p^^. ^ i^p^^.+ie-^f (^'-^i), with Tk = (A;o, e-'^^'k); 

(5) Complex conjugation: (pi. — >■ (.pik; ^p,± — >■ ^-p,±; ^pj ~^-p,3 '^'^^ ^ — > c* , where c 
is generic constant appearing in P{d'^), in V(^) or in {A, J); 

(6. a) Horizontal reflections: (p^ — >■ cri(p'^^^^, (p^ — >■ <^J^k(Ji; ^p,± — >■ ^ii;jp,=F a^c? ^pj — >■ 
-^RhP,rHj(^~'^^^'~^'\ wi^h Rhk = (/Co, -ki, k2) and ml = I, rh2 = 3, r;,3 = 2; 

("^.^j Vertical reflections : (^^ ^ <^Lk; ^p,± ^ ^Rvp,± ".''^d A^j Au^p^rvj, with Ruk = 
(/co, ^1, —^2) o-nd r^l = 1, r^2 = 3, r^3 = 2; 

(^7j Particle-hole : ip^ '^^pk > ~^ '^^Pk> ^p,± ~^ -^-Pp,± (^i^d A^j —A_ppj, with 
Pk= {ko,-ki,-k2); 

(8) Inversion: ^ -^^3^7k,a. ^k,<x ^ -^^lk,<x^3, ^ ^'^S^Zk.a^ 'Pta ^ #lk,a^3, 

^p,± ^ --^ip,± ci'^d Apj ^ ^/p,j, with Ik ^ {-ko, ki, k2). 

Proof. The proof of the fact that P((i\l/), and V(^) arc separately invariant under the 
transformations of the \E' fields has already been discussed in Section 3.1 of [11]. The fact 
that (f)) — (/3L^)~^ o-l^k CT^k (T + ^k CT^k (t) invariant is apparent from the definitions. 
Therefore, here we are left with proving only the invariance of the term (A, J) under the 
transformations (4) to (8) of the hst above. In order to verify these symmetries, it is 
convenient to rewrite the source term in Fourier space: 

i^^J) - E E^P-^*k+p,.^^*k,.+ (B.l) 

ksBt 



^ (^L2)2 ^P,J^k+p,al^'^+^ ^^-^ J^k.CT' 

^ ' peB;3,i i=i,2,3 
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where it is implicit that the terms in the sums with k+p B^^^ should be put equal to zero. 
Symmetry (4)- The term 

(*) := E ^P,rH+P,a^rn,a (B.2) 
k,p,r 

in the first hne of Eq.(B.l) is changed under (4) as: 

k,p,r 

Using the definition of n±, we find that 

Plugging this identity into Eq.(B.3) we see that (*) is invariant under (4). Similarly, the 
term 

(**) := T. ^pM+pA'^+e-'^^''-''^ - icT-e^^''^''^^'^-''^)^^,. , (B.5) 
in the second line of Eq.(B.l) is changed under (4) as: 

(**) ^ J2 ^Tp,,+ie-^^^(^-^^^) • (B.6) 

k,p,i 

Using the definition of (7±, we find that e'^'^+P^^^^-^^^^a±e-'^^^^-^^^ = ^^g±i(k+e)(53-<Ji) ^nd, 
therefore, 

^gi(k+p)(53-5i)^ ^^^_^g-ik(5,-5i) _ .^_^i(k+p)(Sj-Si)y-ik(d3-Si)-^ ^ ^B_y^ 

= e^f (<53-<5i)(i^^e-ik(5,-53) _ ^^_ei(k+p)(5,-53)) _ 

Plugging this identity into Eq.(B.6) and using the fact that k(^j — <53) = (Tk)(<5j+i — di) 
we see that also (**) is invariant under (4). 

Symmetry (5). The term (*) is changed under (5) as: 

(*) ^ E ^-p>-*-(k+p),.^r*:k,. , (B.8) 

k,p,T 

which is the same as (*). Similarly, the term (**) is changed under (5) as: 

(«) ^ 5^(-i-p,)^^(k+p),.( - .a+e+^''(^^-'^^) + za_e-('^+P)(^^-^^))^:,,, , (B.9) 

k,p,J 
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which is the same as (**). 

Symmetry (6. a). The term (*) is changed under (6. a) as: 

k,p,T 

Using the fact that ainT-Ui — we see that this term is invariant under (5). The term 
(**) is changed under (6. a) as: 

(B.ll) 

where 

(B.12) 

Using this identity and the fact that k(6j — ^i) = —{Rh\<:){6r^j — 5i), we see that (**) is 
invariant under (6. a). 

Symmetry (6.h). The term (*) is changed under (6.b) as: 

(*) ^ E ^«.P.-^i(k+p),.^r^R.k,. , (B.13) 
k,p,r 

which is obviously the same as (*). The term (**) is changed under (6.b) as: 

(**) ^ E ^«^P.-^- *i(k+p),.(^'^+e-^"^(^^-^^) - i<7_e+^(»^+P)(^^-^^))^r-^,_^ . (B.14) 

k,p,J 

Using the fact that \<i{6j — 6i) — (i?^k)(6j.^ — ^i) we see that also (**) is invariant under (6.b). 
Symmetry (7). The term (*) is changed under (7) as: 

(*) ^ - E ^-^p.r^;S+p),.^-^Pk,. = E ^-^p.r^Jk,.^r^p(k+p),. , (B.15) 

k,p,r k,p,T 

which is the same as (*). The term (**) is changed under (7) as: 

k,p,i 

= E ^-^P.^- - ^^-e-^"('^-''^))^P(k+p),. ■ (B.16) 

k,p,J 
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Using the fact that k{dj — di) — —{P\c){dr^ — Si) we see that also (**) is invariant under (7). 

Symmetry (8). The term (*) is changed under (8) as: 

(*) ^ E KMo.+^),a hn.(73] *7k,, , (B.17) 

k,p,T 

which is the same as (*). The term (**) is changed under (8) as: 

k,p,i 

Using the fact that aza±az = —cr± and that k(5j — Si) = (/k)(5j — Si) we see that also 
(**) is invariant under (8). ■ 

Appendix C: Symmetry properties of the kernels 

In this appendix we prove Eqs.(2.27)-(2.31). We start by studying the symmetries of 
the kernel quadratic in the external field A and by proving Eqs.(2.27)-(2.29). Next we 
investigate the symmetries of the kernel quadratic in the fermionic fields \1' and linear in the 
external field A and prove Eqs.(2.30)-(2.31). 

The AA kernel. We consider the term quadratic in the external fields in the r.h.s. of 
Eq.(2.25), which has the form (neglecting the dependence on the label M, defining Wtl>(p) 
l^o,2;(tt,b)(p, — p) and assuming, without loss of generality, that l^tt,i,(p) = W\,^f^{—p)): 

E 4,rW^r,r'(P)^-p,r' + 2 5^ ip,rW^rj(p)i-pJ+ A^jWjjip)A_^j> 

P r,T'=+- T=±, jj'=l,2,3 

J=l,2,3 

(C.l) 

which must be invariant under the symmetry transformations listed in Appendix B. Using 
symmetries (4)-(8) of Appendix B, we find that: 

Wr,r'{p) = Wr,ATp)e'^^''-''^^^-^'^ = WIA-P) - 

= W^rM^hP) = Wr,AKp) = WrA^p) , (C.2) 

WrAP) = W-rj+i(Tp)e^§(''-^i)(--^) = -W:/-p) = 

= -W_r,r,j{RhP)e-'''^'^-''^ = Wr,rARvP) = -W^rj(/p) , (C.3) 

WjAp) = Wj+i,j'+i(Tp) = wiA-p) = 

= Wr,,,r,ARhP)e'^^'^-'^'^ = Wr^,,r.ARvP) = W^,y(/p) . (C.4) 
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Eqs.(C.2)-(C.4) imply that one can define natural covariant matrix elements as: 

W,Ap) ■■= e-^'^^'^-'^'^W,Ap) , (C.7) 
which satisfy the following natural transformation rules: 

Wr,r'ip) = WrATp) = W,V(-p) = W^r,-r'iRhP) = Wr,r'iRvP) = W^r,r'(^P) , (CS) 
W^.,,(P) = l^.,,+l(Tp) = -W:A-p) = -W_^,r,ARhP) = W^,rARvP) = -WrAlp) , 
WjAp) - Wj+l,f+liTp) = W*,v(-p) - Wr,j,r,ARHP) = W^w.,'(i?.p) = W^,J'(/P) ■ 

At first order in p, defining Vrjjt,(0) =: ajt, and dp^Wi\,{0) =: fej'^,, with jj, b G {+,—,1,2,3} 
and /i G {0,1,2}, from the first of Eq.(C.8) wc get ar,r' = (ot.t')* = cl-t-t' (i-e., a^-.r' = 
a + a'rr', for some a, a' G M) and h'^,^., = 0, V// G {0, 1, 2}, that is: 

Wr,r'{p) = a + a'rr' + 0{p^) ^ Wr,T'{p) = e-'i^'''^-^'^[a + a'TT'] +0{p^) , a, a' G M , 

(C.9) 

which proves Eq.(2.27). Similarly, from the second of Eq.(C.8), we get a^j = 0, b\.j — for 
Z = 1, 2 and = 6t, for some 6 G M, that is 

W^rj(p)=fcTpo + 0(p') WrAp)-bTpo + 0(p^) , 6GM, (C.IO) 

which proves Eq.(2.28). Finally, from the third of Eq.(C.8), we get that ajj' — cSjj' + c', for 
some c, c' G M and b'^j, = 0, for all // G (0, 1, 2} and j,f G (1, 2, 3}; that is, 

Wj,f{p) = cSj,f + c' + 0(p2) ^ Wr,j{p) = e'^^^^-^^'^ [cSj,f + c'] + 0{p^) , c, c' G R , 

(C.ll) 

which proves Eq.(2.29). 

The Ail^il^ kernel. We consider the term quadratic in and linear in in the r.h.s. 

of Eq.(2.25), which has the form (neglecting the dependence on the label M and defining 
l^„(k,p) :=iy2,i;tt(k + p,k,p)): 

T^E [E4p^£p"A(k,p)C'"+ E 4pC^r>.(k,p)C"] , (C.12) 



43 



which must be invariant under the symmetry transformations hsted in Appendix B. Using 
symmetries (4)-(8) of Appendix B, we find that: 

WriKp) = e^^§(^^-^i)e-*(^+P)(^3-^i)^W^^(rk,rp)e^'^(^3-^^)^ = W;{-k, -p) = (C.13) 
= a,W^r{RhK RhP)(Ti = Wr{R,k, R,p) = W^iP{k + p), -Pp) = a^WrilK Ip)as , 

= -e*P(^^-^^)ail^,,, (i?;,k, P;,p)ai = Wr^R^k, R,p) = (C.14) 
= -Wf (P(k + p), -Pp) = -asWjilk, Ip)as . 

At k = p"^ and p = and defining Wr{pp,0) =: X]j^=o '^^.t^a*) with ao = 1 and Ui, (72,(73 
the standard Pauh matrices, the last identity in Eq.(C.13) reads: 

which imphes a^^ = a^^ = 0; given this fact, the second identity in Eq.(C.13) imphes that 
^°,T — (^-w.t)* ^w,T — (^-w,t)*! third identity imphes that a° is even in r, while 
,^ is odd in r; the fourth and fifth identity imply that both a^, and a^ ,^. are even in uj. 
In conclusion, 

W^r(PF,0) = a° + raV3 , (C.16) 
with a° and two real constants. This proves Eq.(2.30). 

Similarly, defining Wj(pp,0) —: X]^=o ^wj'^m' ^^^^ identity in Eq.(C.14) reads a° j + 
+ (^l,j^2 + aljas = -a° J. + a^^ai + ^jCT2 - aljcr^, which implies that a^j = a^^^- = 0; 
the first identity in Eq.(C.14) reads: 

+ alja2 = e-^'^^^ (ai,,.+i(7i + a^,,.+ia2)e^"^^ (C.17) 

, aiA _• Aim , 
which implies '-^ = e 3 . ^^j^g second identity in Eq.(C.14) implies that 



fl2 / 



a, 



= — ((j^i^j)* and a^j = ((5^?.^^^)*; the third identity implies that a^j = —aj^^r^j 

a^j — o^^ruj'-i fourth identity implies that aj^j = oL,^ j.^^, with / = 1, 2; the fifth identity 

implies that a^j = — ^-wj- Using the second, third, fourth and fifth 

identities for j — 1 immediately gives = and a^ i — a & R. At this point, using the 

first identity, we get 2 = ~'^w,3 — ^0,^ and 2 = ct^.s ~ ~f ' which means 

\/3 1 - \/3 1 

W^i(Pf) = , W^2(Pf) = a{u—ai - -02) , W^3(Pf) = a{-u—ai - -02) , 



(C.18) 



44 



that is 

Wj(p% 0) = ae^'^^(^-^)'^V2 , (C.19) 
which proves Eq.(2.31). ■ 
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